Manifestation of classical bifurcation in the spectrum of the integrable quantum dimer 
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6J[)' We analyze the classical and quantum properties of the integrable dimer problem. The classical 

version exhibits exactly one bifurcation in phase space, which gives birth to permutational symmetry 
■ broken trajectories and a separatrix. The quantum analysis yields all tunneling rates (splittings) 

in leading order of perturbation. In the semiclassical regime the eigenvalue spectrum obtained 
by numerically exact diagonalization allows to conclude about the presence of a separatrix and a 
bifurcation in the corresponding classical model. 
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The problem of correspondence between classical and quantum-mechanical properties of nonlinear systems is cur- 
rently an object of wide interest Q. One interesting topic concerns Hamiltonian systems with a given symmetry (e.g. 
, 1 * ' some permutational symmetry), where classical trajectories exist which are not invariant under the corresponding 
symmetry operation. This topic appears in analyzing selective bond excitation in chemistry and in the quantization 
CZ2 | of discrete breathers [0 . 

We consider an integrable system with two degrees of freedom (TDF), whose classical version exhibits exactly one 
bifurcation (of periodic orbits) and separatrix manifold. This manifold cuts the phase space into three parts - one with 
invariant trajectories, and two with noninvariant trajectories, where the corresponding symmetry is the permutational 
one. By varying a single parameter it is possible to 'switch' between these phase space parts by crossing the separatrix. 
It appears natural to expect in the quantum case a drastic change in the splittings of energy levels (which should be 
zero in the classical limit for the noninvariant phase space parts). However the splittings arc nonzero for any given 

1 value of the control parameter. The only way to avoid contradiction between the classical and quantum cases is to 
assume that the quantum level splittings tend to a step-like function (of e.g. the level pair number) in the classical 

t-H , limit. The step should occur at the position of the classical separatrix. This problem can be also coined dynamical 
tunneling through a separatrix. There exist studies of the influence of classical chaos on dynamical tunneling jj), 
£^ . however we are not aware of any systematic study in the absence of chaos. 

' We are able to trace the splittings of the level pairs using quantum perturbation methods. We consider the 
quasiclassical regime and show that the step indeed occurs. Therefore we are able to extract information about 
the classical separatrix and bifurcation. Further we show, that the quantum density of states (the second integral 
of motion is fixed) exhibits a sharp maximum at the separatrix energy. By calculating the corresponding classical 
quantity (with the help of Weyl's formula) we find that this singularity appears due to the integration over a part of 
the separatrix manifold which includes a hyperbolic isolated orbit. 
Let us consider the integrable dimer model with Hamiltonian 14] 



H = \ {Pf + Pi + X\ + Xl) + 1 ((P2 + X\f + (P 2 2 + Xlf) + ^ (X,X 2 + PiP 2 ) . (1) 



Here Pi, 2, -Xj.,2 are canonically conjugated momenta and positions of two degrees of freedom. System (|l]) is integrable, 
because the classical Poisson bracket of 



B = Pl + P| + Xl + Xf (2) 

un 

With * = l/V2(X + iPY@) becomes 



- 1—1 

X 

with H vanishes. Further (|]J) is invariant under permutation of indices. 



Pr=* 1 * 1 +*^ 2 + i((*^ 1 ) 2 + (^*2) 2 )+C(*t*2+*^l) • (3) 

The equations of motion become ^1,2 = idH/d^\ 2 - 

Isolated periodic orbits (IPO) satisfy the relation grad-ff || gradP. Let us parametrize the phase space of (||) with 
^1.2 = Ai^e 1 ^ 1 ' 2 , Ai^ > 0. It follows Ai, 2 time independent and <f>i = 4>2 + A with A = 0, ir and 01,2 = u> time 
independent. Solving the algebraic equations for the amplitudes of the IPO's we obtain 
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l:A 2 h2 = ±B, A = 0, u = l + C+^B , (4) 
U:A\ = ^B, A=7r, u=l-C+\B , (5) 



III : A\ = -B \1 ± v/l-4C 2 /B 2 J , A = , w = 1 + B . (6) 

IPO III corresponds to two elliptic solutions which break the permutational symmetry. IPO III exist for B > Bb with 
Bb = 2C and occur through a bifurcation from IPO I [Q. The corresponding separatrix manifold is uniquely defined 
by the energy of IPO I at a given value of B > Bb- This manifold separates three regions in phase space - two with 
symmetry broken solutions, each one containing one of the IPO's III, and one with symmetry conserving solutions 
containing the elliptic IPO II. The separatrix manifold itself contains the hyperbolic IPO I. For B < Bb only two 
IPO's exist - IPO I and II, with both of them being of elliptic character. Remarkably there exist no other IPO's, and 
the mentioned bifurcation and separatrix manifold are the only ones present in the classical phase space of (|l|) Jlj . 

To conclude the analysis of the classical part, we calculate the energy properties of the different phase space parts 
separated by the separatrix manifold. First it is straightforward to show that the IPO's (Q)-(^|) correspond to maxima, 
minima or saddle points of the energy in the allowed energy interval for a given value of B, with no other extrema or 
saddle points present p|. It follows 

Ei = H{lYO I) = B + \b 2 + CB , (7) 

E 2 = H (IPO II) = B + jB 2 - CB , (8) 

E 3 = H(lPO III) = B + ^B 2 + C 2 . (9) 

For B < B b we have Ex > E 2 (IPO I - maximum, IPO II - minimum). For B > B b it follows E 3 > Ex > E 2 (IPO 
III - maxima, IPO I - saddle, IPO II - minimum). If B < Bb, then all trajectories are symmetry conserving. If 
B > Bb, then trajectories with energies E\ < E < E3 are symmetry breaking, and trajectories with E 2 < E < E\ are 
symmetry conserving. 

The quantum eigenvalue problem can be properly analyzed in second quantization, which amounts to replacing the 
complex functions 'J, \&* in (||) with the boson annihilation and creation operators a, a + with standard commutation 
relations (to enforce invariance under exchange W ^* the substitution has to be done after rewriting vf"]/* = 
l/2(*4<* + ^**)): 

5 3 1 
H = - + -(afa 1 + a+a 2 ) + -({a+ ai ) 2 + {a+a 2 ) 2 ) +C(a+a 2 + a+ ai ) . (10) 

Note that fi = 1 here, so the eigenvalues b of B — a^ai + a 2 a 2 are integer numbers. Since B commutes with H we 
can diagonalize the Hamiltonian in the basis of eigenfunctions of B. Each value of b spans a subspace of dimension 
(6+1) in the space of eigenfunctions. These eigenfunctions are products of the number states \n > of each degree of 
freedom and can be characterized by a symbol \n, m > where we have n bosons on site 1 and m bosons on site 2. For 
a given value b it follows m = b — n. So we can actually label each state by just one number n: \n, {b — n) >= \n). 
Consequently the eigenvalue problem at fixed b amounts to diagonalizing the matrix H nm with 



| + |fe+i(n 2 + (6-n) 2 ) n = m 

C y/ n(b + 1 — n) n = m + 1 

C- v /(n + l)(6-n) n = m - 1 

else 



(11) 



and n, m = 0, 1, 2, b. Notice that the matrix H nm is a symmetric band matrix. The additional symmetry H nm — 
Htb-n).(b-m) is a consequence of the permutational symmetry of H . 

For C = the matrix H nm is diagonal, with the property that each eigenvalue is doubly degenerated (with exception 
of the state \b/2) for even values of b). The classical phase space contains only symmetry broken trajectories, with 
the exception of IPO II and the separatrix with IPO I (in fact in this limit the separatrix manifold is nothing but 
a resonant torus containing both IPO's I and II). So with the exception of the separatrix manifold, all tori break 
permutational symmetry and come in two groups separated by the separatrix. Then quantizing each group will lead 
to pairs of degenerated eigenvalues - one from each group. There is a clear correspondence to the spectrum of the 
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diagonal (C = 0) matrix H nm . The eigenvalues Hoo = Htb correspond to the quantized IPO's III. With increasing 
n the eigenvalues H nn = H^_ n j^_ n j correspond to quantized tori further away from the IPO III. Finally the states 
with n = b/2 for even 6 or n = (b — l)/2 for odd b are tori most close to the separatrix. Switching the side diagonals 
on by increasing C will lead to a splitting of all pairs of eigenvalues. In the case of small values of b these splittings 
have no correspondence to the classical system properties. However in the limit of large b we enter the semiclassical 
regime, and due to the integrability of the system eigenfunctions should correspond to tori in the classical phase space 
which satisfy the Einstein-Brillouin-Keller quantization rules |Q. 

Increasing C from zero will lead to a splitting AE n of the eigenvalue doublets of C = 0. In other words we find 
pairs of eigenvalues, which are related to each other through the symmetry of their eigenvectors and (for small enough 
C) through the small value of the splitting. Let us calculate the splittings in leading perturbation order. This is 
done by applying standard perturbation theory to each of the states \n) and \{b — n)) and calculating the perturbed 
eigenvectors until the matrix element of the two perturbed eigenvectors with H does not vanish. Due to the band 
structure of our matrix the final result has the following form [|5| : 



AE n = 2 \{ H ii{i+l) J] (H nn -H u ) 1 . (12) 

i=n 

For even b with n = n — b/2 and (O) it follows 



i=n-\-l 



(2|n| - 1)! 2 (| - |n|)i 
For odd b with n = n - b/2 + l/2sgn(ra - b/2) and (111]) we find 



= 2C 2 N — ^±JP — . (13) 



AE n = 2C 2 ^ ( ^ + jgjH 

(2|n|-2)F(*+l-|n|)! 



(14) 



The integer h counts the pairs of equal diagonal elements of ( |ll"|) from the center of H nm towards the corners (b even: 
\n\ = 0, 1,2, ...,6/2 and b odd: \n\ = 1,2, ...,(& + l)/2). Note that for the corner states the obtained expression for 
the splitting is identical with the results in p|. Let us define |n| = ab/2 with < a < 1. For fixed a application of 
Stierling's formula to (|l|),(^4|) yields 



1-aJ 1 ' ' 2a(a6-l) 

For large ab the expression (|l^) should be close to zero if 7 < 1 and its inverse should be close to zero if 7 > 1. So the 
perturbation result predicts a step-like change in the splitting values for 7 = 1 in the limit of large ab. The considered 
asymptotic limit corresponds to the classical limit of ([!(]) . Thus we expect that the splittings of the eigenvalue pairs 
which correspond to symmetry broken classical tori should vanish in this limit. Consequently the condition 7=1 
predicts the position of the classical separatrix with respect to the variable a. 

Now we calculate the eigenvalue spectrum of @ numerically)] (for b = 20 this was done in Q). In Fig.l we 
show the eigenvalues (grouped with respect to their eigenfunctions being symmetric or antisymmetric with respect to 
permutation) as a function of n for b = 600 and C — 50. The classical model has symmetry broken trajectories, and 
a separatrix with energy E sep = E\ = 120600. For the quantum problem we find an inflection point in the eigenvalue 
spectrum of each subgroup at precisely this energy (n w 150). Since n(E) is the integrated density of states, its 
derivative with respect to E gives the density of states p(E), which hereby exhibits a peak at the separatrix energy 
of the classical system (inset in Fig.l). Using Weyl's formula we can calculate its classical counterpart |ij 

Pc i(E,b)= j ' A 2 PA 2 X8{E-H{P,X))8{b-B{P,X)) . (16) 

This integral can be rewritten as p c i(E, b) = § l/dyi/llV-Blsm-O^S, where the integration is done over the surface of 
constant H and B and O is the angle between the two gradients. The denominator vanishes on IPOs. Expanding the 



1 This was done using standard Fortran routines with double precision. When splittings had to be calculated with values 
below 1CP 16 we used Mathematica routines, where the precision can be of any value M. 
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denominator in a Taylor series in the neighbourhood of an IPO it follows, that for elliptic IPOs no singularity develops 
(because the torus surface vanishes) whereas for hyperbolic IPOs (i.e. on the separatrix) a logarithmic singularity 
appears. 

By parametrizing the classical phase space using A\p and <p\,2 the expression ( |l6| ) can be reduced to a single 
integral: 

Pcl (E,b) = - [ dy . (17) 

n J ^C 2 b 2 -W 2 y 2 -{E-b-b 2 /A-y 2 f 

The integration has to be done over all values of y where the expression under the root is nonnegative. This integral 
shows up with a singularity at the classical separatrix energy. The numerical integration is compared in the inset in 
Fig.l with the quantum density of states. We find excellent agreement. In the inset in Fig. 2 the splittings are shown 
with respect to h. The splittings become anomalously small in the region of classical symmetry broken solutions, 
which is bounded again by the separatrix energy. In Fig. 2 we compare the numerically obtained splittings with the 
perturbation theory result (6 = 150, C = 10). Even though the true splittings become as small as 10~ 100 compared to 
the averaged spacings, the perturbation theory reproduces at the best the order of magnitude, but fails by e.g. 50% 
in the absolute value. Consequently we note that higher order terms in the perturbation theory are important even 
when the true splittings are anomalously small. 

Still there is useful information in the perturbation result as shown in ([l5|). In Fig. 3 we show the classical separatrix 
energy E\ for different values of C (b = 600) and compare it to the peak energy in the quantum density of states 
and to the condition 7=1 (which gives us a certain a, which in turn yields a given h and through the numerically 
obtained quantum eigenvalue spectrum a corresponding energy). First we note the remarkable agreement between 
the classical curve and the exact quantum counterpart. But even the perturbation theory gives values which deviate 
by only 6% from the exact result. So while the perturbation theory fails in reproducing the absolute values of the 
splittings, it still contains the information about a classical separatrix with good precision. 

Finally we can easily trace the classical bifurcation by considering the dependence of the largest eigenvalue of the 
quantum spectrum as a function of C: E max = f(C). According to the classical system this function is given by (Q) 
for C > 6/2 and by (|]) for C < 6/2. Differentiating this function twice with respect to C should thus yield a step 
function with the step located at C = 6/2. In the inset of Fig. 3 d 2 //dC 2 is shown for 6 = 600. The step at C = 300 
is nicely observed. 
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FIGURE CAPTIONS 

Fig.l: 

Eigenvalues of the symmetric eigenstates (solid line) and antisymmetric eigenstates (dashed line) versus quantum 
number h for 6 = 600 and C — 50. 

Inset; Density of states for the eigenvalue spectrum from above (solid line) versus energy. The dashed line is the 
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classical prediction using Weyl's formula. 
Fig.2: 

Eigenvalue splittings versus quantum number n for b — 150 and C = 10 (calculated with precision 512). Solid line - 
exact diagonalization, dashed line - perturbation theory result. Note that even for n w 80 the ratio of both values is 
of the order of 0.5. 

Inset: Eigenvalue splittings versus n for b = 600 and C = 50 (compare Fig.l) from exact diagonalization. Splittings 
are of the order of average spacing for n < 150 and collapse to zero for n > 150. 

Fig.3: 

Separatrix energy versus C for b = 600 for the classical system (solid line) . The thick long-dashed line is the position 
of the maximum in the quantum density of states, the thin dashed-dotted line is the perturbation theory prediction 
(7 = 1)- 

Inset: Second derivative of the C-dependence of the maximum eigenvalue of the quantum spectrum for b = 600 versus 
C. The classical prediciton is a step function with values 2,0 and step position C — 300. 
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